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Abstract 

We note that the Thomas Fermi limit of Gross Pitaevskii equation and N >> 1 limit of quantum mechanics, where N is 
the dimensionality of space, are based on the same point of view. We combine these two to produce a modified Thomas 
Fermi approximation which gives a very good account of the energy of the condensate in harmonic trap. 
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I. Introduction 

Bose Einstein condensation has been experimentally 
achieved with particles in a trap-generally a simple har¬ 
monic oscillator potential. The condensate is well de¬ 
scribed in terms of the Gross Pitaevskii equation (GPE)P, 
iii both for the energy (chemical potential) of the sta¬ 
tionary state and as well as dynamics. To find the lowest 
energy of the system, one can use a variational estimate 
with a Gaussian trial function for a Thomas Fermi ap¬ 
proximation which ignores the kinetic energy of the par¬ 
ticles. The Thomas Fermi approximation works well for 
a large number of particles but generally is an underes¬ 
timate as compared to the variational calculation. Hence 
finding the correction to the Thomas Fermi limit is of inter¬ 
est. The first attempt at doing this is the work of Schuch 
and Vin who combined a WKB type of analysis with 
the Thomas Fermi idea to obtain a very effective modified 
Thomas Fermi energy. It has not been noticed yet that 
analogous to the Thomas Fermi limit is the infinite dimen¬ 
sion (large N) limit for finding th^round state energy of a 
quantum mechanical problem [a, IT], [sj, [Oj. As N ^ oo, the 
kinetic energy becomes negligible and one finds the ground 
state energy from a minimization of an effective potential. 
Corrections are obtained in a systematic fashion by first 
considering simple harmonic motion about the minimum 
and then including anharmonicity. Our point in this arti¬ 
cle is that the first correction to the Thomas Fermi energy 
can be found in a manner analogous to the correction to 
the A —>■ 00 limit and also the large N approach can be 
combined with the Thomas Fermi limit to yield accurate 
answers more easily. 

In sec II, we will recall the large N quantum mechani¬ 
cal problem. We will first consider the harmonic oscillator 
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potential where the exact answer is obtained at the leading 
order itself and all subsequent corrections are zero. More 
pertinent is the an-harmonic oscillator which we treat next 
and show that a two term answer is a significant improve¬ 
ment on the leading terms. In section III, we use the phi¬ 
losophy of approach to find the correction to the Thomas 
Fermi limit in a one dimensional problem. In sec IV, we 
combine the large N approach with the Thomas Fermi 
limit to obtain the large N approach with the Thomas 
Fermi limit to obtain a very reasonable estimate of the 
stationary state energy. 


II. 1/N expansion 

In this section, we recall the 1/A expansion for the har¬ 
monic oscillator and the an-harmonic oscillator. We will 
get the ground state energy to two term accurately in each 
case and compare with the known exact (analytic or nu¬ 
merical) answer at A = 1, which is the worst situation for 
the technique.Even in this very unfavorable situation, the 
comparison between approximation and reality is good, 
giving us confidence in this procedure. 


ILL Harmoriic oscillator potential 

The laplacian operator for the zero angular momentum 
state in a N-dimensional space has the form 




^ (j.N-1 ^ \ 

r^-1 dr^ dr 


( 1 ) 


Consequently, the Schrddinger equation for the ground 
state of the simple harmonic oscillator of frequency w reads 


h? \ d 
2m r^~^ dr 


N-l 


(r^-'T^) + ^mccV^ 
dr 2 


Ecj) 


( 2 ) 


where, m is the mass of the particle and E is the ground 
state energy. If we carry out the transformation (j) = , 

r~^ 
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then 


Considering the scaling as 


Hu=- 




2m ^ 
= Eu 


(TV - l)(iV-3) ,1 2 2 

--- -u] H— muj'^mu 

4r2 2 


( 3 ) 


where, ^H' is the Hamiltonian operator. We define a di¬ 
mensionless spatial variable ^ as 


m\ 

E = 7V4/3(^)2/3;^1/3 

m 

Eq. ® reduces to 


( 10 ) 

( 11 ) 


= rV- N 

muj 


so that the Hamiltonian now turns out to be 


( 4 ) 


H = N\-^^ + (1-^)(1-^) , 
feu ^ 2m de 8^2 2 


( 5 ) 


Rescalling E = NHuie, we can rewrite Eq. ® 

_+ (l-^)(l-f) ^.2 E 

2m 8^2 +2^ iVfeu“ (6) 

= eu 


As —>■ oo, only the 0(1) terms on the right hand side of 
Eq. m survive and in the ground state the particle sets 

1 

at the minimum of the effective potential V = ^ -I- 
The minimum is at ^ with = 1/2 and in this 

approximation the ground state energy is eo = 1/2 and 
given by 


Eq = Nfiuieo = 


Nfvjj 

2 


( 7 ) 


This is what in critical phenomena would be called the 
spherical limit Eliil- In this context it is the leading 
term for >> 1. The first correction to the leading 
order is obtained by expanding V about and writing 
iJ = ^ + E', where ^ is the N ^ oo form and for finite 
N, El' is the variation around the infinite N limit. To the 
leading order, H' can be expressed by expanding around 
^0 = 1/2 


H' = 


1 d^ 

2m dc 


+ 2(?-Co)" 


1 

N 


( 8 ) 


which is the Hamiltonian of simple harmonic oscillator of 
frequency oj' with cu'^ = i/d^ and an additional part. The 
ground state energy is clearly ^ = 0 and hence there 

is no correction at this order. The ground state energy 
thus remains Eq = + 0(;^) which is comparable to 

the exact answer of Higher order corrections will 

remain zero because in this case is the exact answer. 


Hu = - 


1 d^u 


(1-^)(1-^) 1 , 

--"4] + —u = eu (12) 


In the limit of TV —>■ oo, there is only the effective potential 
E = ^ -I- ^ and for the ground state, the particle sits 
in the minimum of the potential which occurs at ^ = ^O) 
where = 4“^ and Enm = ^ x 4^^^ giving eg = ^ x 4^/^. 

To obtain the first correction H', we need to keep the 
derivative term in the Hamiltonian of Eq. m- Keeping 
the term — „ E^ in Hq and expanding the effective poten- 
tial V about ^g to quadratic order we get 


E~ 4 X 4^/3 2.45(A^)2 


In this approximation 


H' = 


I 

d(A^)2 


-h 2.45(AO^ 


41/3 

w 


The ground state energy ei of TV' is 

_ 1 4I/3 ^ 0.31 

~ 2 y m~'m ~ ~ 

To this order the energy E is given by 
E = )V=a1/3„4/3|^ ^ 4>/» + /II + 0(4^)1 (13) 


For A = 1, at this order we get E ~ 0.61('^)^/^A^/^, 
which is to be compared with the numerically exact value 
of E ~ 0.668(^)2/3 a 1/3^ Considering the fact that A = 1 
as far remote form of A 00 as possible and we have 
calculated to a reasonably trivial order, this agreement is 
quite impressive. The agreement of the two term answer 
above is far more impressive for A = 3. In the next section, 
we consider the Gross Pitaevskii model 


III. Gross Pitaevskii Model with a harmonic trap 
in IN 


11.2. Anharmonic oscillator (Xx'^ potential) 

Had the potential been Ax"^ in place of ^muPx^, we 
would have been able to formulate an identical scheme of 
expansion and corrections could no longer be zero. In this 
case, Eq. (fS]) would take the form 


Hu=- 




2m ^ dr'^ 
= Eu 


-5x-“I+ r “ (9, 


For the N particle Bose Einstein condensate in a har¬ 
monic trap in IN, the wave function ‘u’ satisfies the Gross- 
Pitaevskii equation (GPE) 


ll d U 1 22 ||2 7-? 

— -- — + -mu! X u + g\u\ u = Eu 

2m dx‘‘ 2 


(14) 


where, the coupling constant ^g' corresponds to the two 
body interaction between the particles that constitutes the 
condensate. For g > 0, it signifies repulsive interaction 
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while for 5 < 0, attractive interaction is indicated. In this 
work we will be concerned with g > 0. We carry out the 
following rescalings, 


x= \l -^ 

muj 


(15) 


q = Xhjjj\ - 

V muj 

to arrive at the following 


1 (Pfj) 

2de 


2 




The scaling given in Eq. (HU), ensures that 


(16) 


= 1 . 


III.l. Energy obtained from variational calculation: 

Considering the trial wave function of the form, (j){x) = 
/ 2 A energy of the IN condensate is ex¬ 
pressed as 


^ 1/1 a2n 1 


(17) 


Minimizing with respect to A, we get the value Aq at the 
minimum to satisfy 


1 


2A3 


-b 


A 


\/^Aq 


Aq 

2 


(18) 


Obtaining Aq for different values of A, we can find the 
ground state energy from Eq. dm as a function of A and 
this is plotted as the solid curve in Fig. [5] 

III.2. Energy obtained from large A expansion 

To explore the large A behavior we consider the scale 
transformation ^ = A^/^y 


1 d?(b , , ,,2 , 


(19) 


where E = eXfuo. We would like to compare the above 
with Eq. The present situation is more complicated 
since instead of a term proportional to we now have 
which can be known only after the problem is solved. 
However we note that |(()p —?> 0 for ^ ^ 1 and hence the 
two terms are of the same character. We try to exploit 
this analogy. 

We take the limit of A —>■ 00 and drop the first term in 
Eq. m to write the solution (p = 4>o{y) at this order as 


\ 4 >o\^ = Co — 
= 0 


for y <yl, = 2 eo 
otherwise 


( 20 ) 


This is the Thomas Fermi approximation and the value of 
ep is fixed by the normalization condition = f |<(>|Qdy = 
1 ( 260 )^^^^ which leads to 


1 / 3 
2 4ai/2 


)2/3 


( 21 ) 


If we view the Hamiltonian as H = — + \ + , 

then we have taken into account the part o^ H which is 


\ -\-\(p\’^ and this extends over the region |y| < \f2^. The 

1 ri^ 

remainder of the Hamiltonian is — 

2A^ dy^ 


over all space 

and a potential which is ^ for |y| > v' 2 eo and ep for 
I 2 /I < V2eo- The lowest energy eigen value of the potential 
is ei relative to a base value of ep. and the corresponding 
V(y) is shown in Fig. [TJ Hence we need to solve 


1 d^cj) 

~ ^d^ 

With the analytic form of V{y), as 


fo(y)0 = (eo + ei)0 


y{y) = 


r 

2 

= eo 


for |y| > ^/^ 
for |y| < 


( 22 ) 


(23) 


ViV) 



Figure 1: Above figure describes the schematic diagram for the po¬ 
tential given in Eq. ll23ll . All other parameters that we have con¬ 
sidered in the further analysis of WKB have also been shown in this 
figure. 

It will be very cumbersome to find the exact answer 
(which can be done by matching trigonometric functions 
and Weber functions at y = ±^/2^) and so we resort to a 
WKB procedure which for the lowest value of ei yields 

2A[ f \/2(ei -b ep — V{y))dy+ 

° (24) 


/o 

rvo 


2 (ei + ep - Y)rfy] = I 


Where yp is to be obtained from y^ = 2(ep -b ei). The first 
integral on the left hand side of Eq. [M] is easily seen to 


3 



















be 2^eoei, while the second integral is evaluated to the 
leading order in ,/—. We have 


[ \/2(ei + eo) - y'^dy = (ei + eo)[6»o - ] (25) 

Jo Co + Cl 


where, Cos^Oq = Extracting Oq to we get 


/n = =(^-1)^ 

Co 4 Aeo 

“ 2^^ A2/3 


(26) 


We have the two term Thomas Fermi energy consequently 
given by 




(27) 



with r scaled by 



(the oscillator length) and we get 


1 N muj. N --1 g \u\^u 


E 

lu^' 


(30) 


with the normalization condition f \u\'^d^ = ■ We 

now make a further rescalling by N as done in sec II and 
with ^ we reach into the following 


I (Pu 
2iV2 dR^ 


u R^ , moj , N-i g 

8r5 + -“+(—) ■ 

luP E 


(31) 


rN-1' 


hjjjN 


In the large g, large N limit keeping the ratio 
finite, the term O(^) can be dropped and we have the 
modified Thomas Fermi form for the wave function given 

by 


wo; g |np 

^ fi ’ 




N-1 


(32) 


where e = and the normalization now becomes 


(33) 

F(iV/2) V mw 

where Ri and R 2 are the limits between which u is non 
zero i.e., i?i and i ?2 correspond to the zeroes of the right 
hand side of Eq. with 


Rl,i = e ± ^^2 - - (34) 

The normalization condition given in Eq. (1331) becomes 
with the wave function of Eq. (15^ 


Figure 2: Black curve shows the ground state energy of IN GPE, 
obtained after including the first order correction over the usual 
Thomas Fermi approximation. Variational result is indicated by the 
red curve and the blue curve shows the energy values obtained from 
TF approximation. 


27r^/2 


^ ^ '^N/2 


r(iv/2) TOW 

(i?^+2 - i?f+2) 
2(7V + 2) 


Vf+1 r e 

9 Lv 

{R^-^ 


{R2 -Ri)- 


-R^-^y 


(35) 


8(IV- 2) 


= 1 


IV. Combining large N and Thomas Fermi approx¬ 
imation 


In this section, we seek an improvement on the Thomas 
Fermi technique by starting with the GPE in N dimension 
and doing simultaneously a large N and large coupling 
constant approximation. To this end we begin with 

P 1 

-H —muPPy + giy^^y = Ey (28) 

2 to 2 

as the stationary state GPE and use the transformation 
y = r(jv-i )/2 to write for N » 1 as 


fi2 ^2 
2m dr2 


7V2 1 „ „ |w|2?i 

^]u+-TOW r u + 5 ^ 


Eu (29) 


Along with the expression for Ri and i ?2 from Eq. dSl , we 
have the large N modified Thomas Fermi approximation 
for E in Eq. ((551) . If we now set N=3, after few simple 
algebra we reach into 




Ll5^Sa;^ ^hio'5 10 

9 


(36) 


^ ^ 3/2 

Airhuj h 


If a is the scattering length and n be the total number of 
particles then = :^ an and writing = Qosc, we 

can rewrite Eq. (1361) as (with E = ^) 


4(2i5-3)‘/^[fi» + 5fi-fl 


na 

^osc 


(37) 
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We have found E coming from this equation for 
atoms for which a = 3 x 10“^ agrees reasonably well to the 
value of the ground state energy obtained from standard 
quantum mechanics calculation and the value at TF Fermi 
limit. We have shown in Table 1 the change of ground 
state energy (^) with the total number of atoms (n), 
1 /N correction method having higher accuracy for higher 
values of n. 



n 

Etf+1/n 

Eqm 

Etf(M odified TF of [a) 

200 

1.774 

1.688 

1.642 

600 

2.046 

1.927 

1.877 

1000 

2.245 

2.134 

2.071 

2000 

2.62 

2.535 

2.457 

4000 

3.159 

3.112 

3.025 

6000 

3.571 

3.550 

3.461 

8000 

3.914 

3.914 

3.825 

10000 

4.214 

4.231 

4.142 

12000 

4.483 

4.513 

4.426 

14000 

4.727 

4.770 

4.684 

16000 

4.954 

5.007 

4.921 

18000 

5.165 

5.228 

5.143 

20000 

5.363 

5.435 

5.350 


Table 1: Ground state energy is calculated through l/N expansion 
method in the second column.The frequency of the harmonic oscil¬ 
lator is iu = 207rs“^ and mass of Cs is taken 133 amu. The energy 
values given in 3rd (from general quantum mechanical calculation) 
and 4 th column are taken from Table II of the work by P. SCHUCK 
and X. VINAS @] to compare the accuracy of our results. With in¬ 
crease of n, l/N method gives quite impressive result than TF limit 
as expected. 


Figure 3: This figure has been taken from the book ‘Bose-Einstein 
Condensation in Dilute Gases’ by C. J. Pethick and H. Smith. (Cam¬ 
bridge university press) chap 6 page no. 155. The figure shows vari¬ 
ational estimate of the energy per particle for an isotropic harmonic 
trap as a function of the dimensionless parameter Na/aosc, where ’a’ 
is scattering length, aosc is oscillator length and for this figure only 
N indicates the total number of particle i.e., N = n. The dotted line 
is the result obtained from the Thomas Fermi approximation. 


In general the Thomas Fermi answer for the energy 
of the condensate falls significantly below the variational 
calculation (generally close to the real answer) for small 
values of n because of the neglect of kinetic energy term. 
This is seen from Fig. [3] taken from the work by Pethick 
and Smith Q. 
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Figure 4: Above figure is plotted taking the values (2nd, 3rd and 
4th column) of Table 1 with . This figure clearly signifies the 
strength of large-N quantum mechanics applied in higher dimension 
also. 


V. Conclusion 

In conclusion, we have proposed the method of large-N 
quantum mechanics and have applied this method to vari¬ 
ous ID systems (harmonic oscillator, an-harmonic oscilla¬ 
tor, GP model for dilute BEG at T = 0 both in lower and 
higher dimension) and have derived the significant correc¬ 
tion in the leading order. For all these cases, the correc¬ 
tions have shown sufficient improvement over the base val¬ 
ues. In case of GP model, the energy corrections obtained 
by this method have shown quite a remarkable improve¬ 
ment over the usual Thomas Fermi approximation. Gom- 
parison with the numerical values obtained for the ground 
state energy as given in column 2 of Table 1 justifies the 
potential of this method. 
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